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Abstract: In this paper we consider (n + l)-dimensional cosmological model with scalar 
field and antisymmetric (p + 2)-form. Using an electric composite Sp-brane ansatz the 
field equations for the original system reduce to the equations for a Toda-like system with 
n(n — l)/2 quadratic constraints on the charge densities. For certain odd dimensions 
(D = 4m + 1 = 5,9,13,...) and (p + 2)-forms (p = 2m — 1 = 1,3,5,...) these algebraic 
constraints can be satisfied with the maximal number of charged branes (i.e. all the 
branes have non-zero charge densities). These solutions are characterized by self-dual or 
anti-self-dual charge density forms Q (of rank 2m). For these algebraic solutions with the 
particular D, p, Q and non-exceptional dilatonic coupling constant A we obtain general 
cosmological solutions to the field equations and some properties of these solutions are 
highlighted (e.g. Kasner-like behavior, the existence of attractor solutions). We prove the 
absence of maximal configurations for p = 1 and even D (e.g. for D = 10 supergravity 
models and those of superstring origin). 
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1. Introduction 

In this paper we investigate composite electric S'-brane solutions (space-like analogues of 
D-branes) in an arbitrary number of dimensions D with scalar field and (p + 2)-form. The 
(p+2)-form is considered using a composite electric ansatz and the metric is taken as diagonal. 
All ansatz functions for the metric, form field and scalar field are taken to depend on only one 
distinguished coordinate which is taken as time-like for the cosmological solutions considered 
in this paper. Previously, related work on cosmological and S-brane solutions can be found 
in [1]-[11] and references therein. 

The procedure we use to investigate our system of multi-dimensional gravity plus scalar 
field plus form field is similar to the approach used in [12]. This work also studied a system 
with scalar fields and antisymmetric forms field defined on the manifold Mq x Mi x . . . M n 
(Mi are Einstein spaces and i > 1). The form fields were taken in the form of a composite 
electro-magnetic p-brane ansatz, the metric was block-diagonal, and all scale factors and fields 
depended upon coordinates of Mq. Under these conditions the original model could be reduced 
to a gravitating, self-interacting sigma-model on Mq with quadratic "constraints" on the 
charge densities. These constraints came from the non-diagonal part of the Einstein-Hilbert 
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equations. It was shown that the constraints could be satisfied for certain "non-dangerous" 
intersection rules of the branes [12]. In the present work we use the same sigma-model 
approach to show that it is possible to satisfy the constraints maximally (i.e. all the branes 
carry non-zero charge densities) in certain odd dimensions. We then examine cosmological 
solutions in these odd dimensional cases, and discuss some of their interesting features such 
as their Kasner-like behavior. 

The importance of studying solutions with "maximal" number of branes is related to 
research of oscillating behavior of cosmological solutions near the singularity [13, 14, 15], e.g. 
using the so-called billiard approach [13]. In [14] and other related works (for a review, see 
[15]) it was argued that in superstring cosmology one gets chaotic behavior in terms of the 
"oscillations" of Kasner parameters as one approaches the cosmological singularity. In these 
works the maximal number of electric branes were considered. 

In the next two sections we will give the set up for the system of D = n + 1 dimensional 
gravity, plus a scalar field, plus a (p+2)-form field. For the conditions considered in this paper 
(diagonal metric, composite Sp-brane ansatz for the antisymmetric (p + 2)-form field, and all 
the ansatz functions depending only on one coordinate) this complex system can be reduced 
to a 1-dimensional a-model. This transformation greatly helps in studying the solutions of 
the system. 

In section 4 we consider the quadratic constraints for the charge densities of the branes. 
We find that these constraints have "maximal" solutions with all non-zero brane charge 
densities in particular odd dimensions with particular form fields: D = 5, 9, 13, . . . and p = 

1, 3, 5 . . ., respectively. We prove also the absence of maximal configurations for p = 1 and 
even D. 

In section 5 we investigate cosmological solutions to the field equations for these odd 
dimensions. We look at the proper time behavior of the simplest of these solutions. We also 
show that certain solutions exhibit Kasner-like behavior at these early times. 

2. D-dimensional gravity coupled to g-form and scalar field 

Here as in [12] we consider the model governed by the action 



where g = gMNdz <8> dz N is the metric, ip is a scalar field, A £ 1 is a constant dilatonic 
coupling and 



S= f d D z^/\g~\ R[g}-g MN dM^d N ip-^exp(2X^)F 2 , 




(2.1) 



F = dA= -F Ml M Jz Mi A ... A dz M \ 
n\ q 



(2.2) 



is a (/-form, q = p + 2 > 1, on a D-dimensional manifold M. 
In (2.1) we denote \g\ = \ det(gMN)\, an d 



FM 1 ...M q FN 1 ...N q 9 



,MiJVi 



...g 



,M q N t 



(2.3) 
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The equations of motion corresponding to (2.1) are 



Rmn — tt.'/.w.y R = Tmn, 



A 



,2\<p F 2 = 0j 



V Ml [g](e 2X ^F M ^) = 0. 



(2.4) 

(2.5) 
(2.6) 



In (2.5) and (2.6), A [g] and \/[g] are Laplace-Beltrami and covariant derivative operators 
corresponding to g. Equations (2.4), (2.5) and (2.6) are, respectively, the multidimensional 
Einstein-Hilbert equations, the "Klein-Fock-Gordon" equation for the scalar field and the 
"Maxwell" equations for the (/-form. 

The source terms in (2.4) can be split up as 



Tmn = T MN [<p, g] + e 2Xlp T M N [F, g] , 



(2.7) 



with 



r MN[F,g] = 



Tmn[<P,9] = d M <pd N ip - -gMNdp(pd P (p, 
1 



1 772 , 77 77 M 2 ...M q 

-gUN? + q^MM 2 ...M q ^ N 



being the stress-energy tensor of the scalar and (/-form, respectively. 
Let us consider the manifold 



(2.8) 
(2.9) 



M = (u-,u + ) x R n (2.10) 
with the metric taken to be diagonal and of the form 

n 

g = we *rMdu ®du + Y^ eWMetdy* dy\ (2.11) 

i=\ 

where w = ±1, and u is a distinguished coordinate. The metric ansatz functions j(u), (j) l {u), 
the scalar field <p(u) and the g-forms are assumed to depend only on u. For definiteness one 
can think of u as the "time" coordinate but when w = +1, u is space-like. Here 

Si = ±1 (2.12) 

are signature parameters with i = 1, . . . , n. When u is time- like and all £j = 1 the solutions 
are cosmological. The functions 7, (f) 1 : (u-,u+) — > M are smooth. 

In order to deal in a general way with the different possible indices for the various forms 
((/-form, volume form) we define 

n = {0, {1}, . . . , {n}, {1, 2}, . . . , {1, . . . , n}} (2.13) 
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which is the set of all subsets of 

Io = {l,.-.,n}. (2.14) 

These sets indicate the number and ranges of the indices of the q- forms. 

For any I = . . . , i^} £ 0,q with i\ < . . . < i^, we define a form of rank d{I) = k 

t(I) = dy h A ... Ady ik , (2.15) 

The corresponding brane submanifold is described by coordinates y n , . . . , y %k . We also define 
the £-symbol as 

£{I) = e il ...e ik . (2.16) 
We adopt the following electric composite Sp-brane ansatz for the field of the (p+2)-form 

F=Y^ ^Ar(/) (2.17) 

where the set 

n e = {I €Qo\d(I) = q-l=p+l} (2.18) 

contains all subsets of S7o of the "length" p+1, i.e. of the form {io, i±, i p }. 

We assume that the scalar potential and the scalar field only depend on the distinguished 
coordinate 

& = ^(u) , <p = <p(u). (2.19) 

3. (T-model representation with constraints 
3.1 cj-model 

The system of the previous section can be greatly simplified. In [12] (see Proposition 2 in [12]) 
it was shown that the diagonal part of Einstein equations (2.4) and the equations of motion 
(2.5)-(2.6), for the ansatz given in (2.11), (2.17)-(2.19), are equivalent to the equations of 
motion for a 1-dimensional a- model with the action (see also [2, 3]) 



--hi 



duM 



GijW +V 2 +Y1 £ ( I ) exp[-2f/ / (<A, v?)]^ 7 )' 
ien e 



(3.1) 



the overdots represent differentiation with respect to the distinguished coordinate, i.e. -jj^. 
The factor N is the lapse function given by 

AT = exp( 7o - 7) > (3.2) 

with the definition 

n 

i=i 
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Next the factor in the exponent is given by 

U 1 = U I (cf>,<p) = -\ip + J2<t> i - (3-4) 

Finally, 

Gij = Sij - 1 (3.5) 

are components of the "pure cosmological" minisupermetric matrix, i,j = 1, . . . ,n [16, 17]. 
In this rewriting of the system the generalized "Maxwell equations" of (2.6) become 

4- fexp(-2f/ / )6 /N ) = 0. (3.6) 

They can be readily integrated to give 

¥ = Q(I)exp(2U I ), (3.7) 

where Q(I) are constant charge densities and I E Q e . 

We will analyze the c-model representation of the original system in the harmonic gauge 
where 

7 = 70, AA = 1. (3.8) 

We now further simplify the a- model in (3.1) by introducing collective variables x = 
(x A ) = ((f> l ,(p) and a "truncated" target space metric 

G = G A Bdx A <g> dx B = Gijdft ® dcfP + d<p <g> dip, (3.9) 

{0ab) = (°o °i) ■ (3 - 10) 

U J (x) is defined in (3.4). It can be written as U I {x) = U^x 

(U I A ) = (5 iI ,-X) , (3.11) 

where 

is an indicator of i belonging to /. For fixed charge densities Q(I), I G fi e , the equations of 
motion for the a- model in (3.1) are now equivalent to the equations from the Lagrangian 

Lq = X -Gabx A x b - Vq, (3.13) 

with the zero-energy constraint 

Eq = )pABX A x B + V Q = 0. (3.14) 

Here 

V Q = \Y1 exp[2C/ 7 (x)]. (3.15) 

In section 5 we will examine explicit solutions of the field equations that result from 
(3.13) - (3.15). 
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3.2 Constraints 

Due to diagonality of the Ricci-tensor for the metric (2.11) the non-diagonal part of the 
Einstein equations (2.4) reads as follows 

Tij = 0, i+ j. (3.16) 

This leads to constraints among the charge densities Q(I). First, the non-diagonal components 
of stress-energy tensor are proportional to 

e 2 ^F lM2 ... Mq F^- M \ (3.17) 

with i / j. From (2.17), (3.7) and (3.4) we obtain for the (p + 2)-form 

F = ( p+1 y Qio...i P exp(2^° + . . . + 2</> i * - 2\^)du A dy io A ... A dy** (3.18) 

Inserting this in (3.17) we are led to the following constraint equations on charge densities 
[12] 

n 

Cij = ^ ^ Qiii...i p Qjii...ip^ii ( i ^ • • • £j p 6 ^ =0, (3.19) 
ii : ... : ip=l 

where i ^ j; i, j = 1, . . . , n. is proportional to exp(— 2\ip — 2j + 2(j) 1 + 2<pi)Cij for i / j. 
Herep= q — 2 and Qi^.Ap are components of the antisymmetric form ofrankp+1 = q— 1 

and 

Qioh...i p = Q{{io,h, • • • , ip}) (3.20) 

for iq < i\ < ... < i p and {io, ii, . . . , i p } £ r2 e . The number of constraints in (3.19) is 
n(n — l)/2. In the next section we will show that these constraints can be satisfied when the 
dimension of space-time takes certain odd values. 

4. Solution to constraints on charge densities in various dimensions 

The constraints (3.19) can be rewritten as 

n 

i\ ,... : ip=l 

i 3'i hJ = 1; • • • > n - The charge densities have been redefined via 

v 

Qi h...i p = Qi 0il ...i p n exp(0 ifc ) (4.2) 

k=0 

and the indices were lifted by the flat metric 

r, = eidy 1 ® dy 1 + . . . + e n dy n dy n = r, ab dy a ® dy b , (4.3) 
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where (77^5) = (i] ab ) = diag(ei, . . . ,e n ), i.e. Q 1 °-' 1 p = rf^'o . . . rf^pQ^ . ( Here Cf = 
C{ exp(^ + <p) with C{ = C ik r] k i.) 

These Qion-i P can be viewed as "running" charge densities with the functional depen- 
dence coming from (f) %k (u). The charge densities will vary with time or spatially depending 
on whether u is a time-like or space-like coordinate respectively. 

4.1 Maximal configurations for dimensions D = Am + 1 

D=5 case: To help illustrate the preceding general analysis of the constraints in (4.1) we 

consider the explicit example D = 5, n = 4 and £\ = . . . = £4 = 1. The constraints of eqs. 
(4.1) read 

Q13O23 + O14O24 = 0, (4.4) 

Q12O23 - O14Q34 = 0, (4.5) 

Q12O24 + O13Q34 = 0, (4.6) 

Q12Q13 + Q24Q34 = 0, (4.7) 

Q12O14 - O23Q34 = 0, (4.8) 

QisQu + Q23Q24 = 0. (4.9) 

It is not difficult to verify that the only non-zero solutions to eqs. (4.4)-(4.9) are 

Q12 = TQ34, (4.10) 

0l3 = ±<524, (4.11) 

Qu = TQ23- (4.12) 

This may be obtained by considering the following three pairs of equations: (i) (4.4) and 
(4.9); (ii) (4.5) and (4.8); (iii) (4.6) and (4.7). The solution (4.10)-(4.12) may be written in a 
compact form as 

Qioh = i^' ' 1 ^^ 1 ^ '" = = '-(*^5)*on' (4-13) 

where * = *[r]] is the Hodge operator with respect to ij. That means that any self-dual or 
anti-self-dual 2-form is the solution to a set of quadratic equations (4.4)-(4.9). 

D= 4m +1 case: We now look at the general case. Based on the D = 5 case we will 
take the "running" charge density form Qi ...i p as self-dual or anti-self-dual in order to satisfy 
the constraint equations (4.1). This form can only be self-dual or anti-self-dual when the 
number of the non-distinguished coordinates is twice the rank of the form: 

n = 2(p+l). (4.14) 

Thus, we restrict ourselves to the case when 



Qio...i P = ± (^~7)T £ *o • ■ ■%> Jo ■ ■ -ip Q j °" ' jp = H*Q)io-..i P - (4-15) 
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Here the symbol * = *[rj\ is the Hodge operator with respect to 77. Squaring the Hodge 
operator gives 

(*) 2 =£(-1)(p +1 ) 2 1, (4.16) 

where 

£ = e 1 ...e n (4.17) 

equals ±1 depending on if there are an even or odd number of time-like coordinates. 

It can be easily verified that the set of linear equations (4.15) has a non-zero solution if 
and only if 

£(_1)(p+i) 2 = 1. (4.18) 

The dimension of the space of solutions is ^C^p+i)- ^he factor of \ comes from (anti-)self- 
duality condition. 

We will now demonstrate how having self dual or anti-self dual charge density form results 
in the constraints (3.19) being satisfied. First, consider 

n 

Cj 3 = ^ ] Qiii...i p Q J 1 P = 
U,...,ip=l 

n n 1 

= £ E ±7- T W^ 1 -..Wo.., P JO --^0^ 1 --^, (4.19) 
u,...,i p =i j ,...,jp=i 

where i 7^ j, and we have used the requirement of self duality or anti-self duality for the 
charge density. This can be further rewritten as 

n n 1 
C j = \ \ " ± — r Q33l-3pQ3h-ip 

n,...,ip=l ji,...j p =l 

n n (— -\\P 

\" \ " -L ^ J_ p 03il-ipr>33\-3v 

~ Z_> ' ■> p\ l ^-3p3^i-tp^ V 

i 1 ,...,ip=lj 1 ,...,jp=l 

= {-VPC*. (4.20) 

Note that j is not summed over in the two sums above. In going from the second line of 
(4.19) to the first line of (4.20) we have carried out p+1 identical sums with: jo = j, ji = j, 
j p =j, respectively. 

From (4.20) one finds that the constraints in (4.1) are satisfied automatically for odd 
p = 2m — 1 since 

<v = -<v =► <v = . 

From (4.18) we get for odd p 

6 = 1, (4.21) 
i.e. the metric (4.3) is either Euclidean or has an even number of time-like directions. 
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The relationship between n and p is n = 2(p + 1). Thus, we find non-trivial solutions to 
the constraints (4.1) when the total spacetime dimension is 

D = n + l = 4m + 1 = 5,9,13,... (4.22) 

and the signature parameter £ is positive. 

4.2 Absence of maximal configurations for p = 1 and even D 

In this subsection we show that for the even dimensional case (D = 2k) with 3-form (p = 1) 
there are no solutions with maximal number of electric 51-branes. 
Here we put e± = . . . = e n = 1. Equations (4.1) imply in this case 

n 

^ ^ QiiiQjii = $ijPi = Pij- (4.23) 
h=l 

The indices are not summed in the second term. Now we assume that all Qij / 0, i / j, 
i.e. a composite 1-brane configuration with maximal number of electric branes ("strings") 
is considered, and show that this leads to an inconsistency. The constants Pi, which are the 
values of an n x n diagonal matrix, satisfy Pi > 0, since when i = j the first equation is just 
a sum of squares. The exact values of Pj's will not be needed here. In matrix notation (4.23) 
reads 

-Q 2 = P, (4.24) 

where we have used the antisymmetry relation Qij = —Qji- Calculation of the determinants 
of the matrices in the previous relation leads to 

(-l) n (detQ) 2 = detP > 0. (4.25) 

which is not valid for odd n. (For odd n detQ = 0.) Thus, there are no "maximal" solutions 
to constraints for even dimensions D and p = 1 in the model under consideration. 

This implies the absence of "maximal" configurations of composite electric 51-branes 
in 10-dimensional supergravities and low-energy models of superstring origin when only one 
3-form is considered. 

In the next section we examine the cosmological type solutions to the field equations for 
D = 5, 9, 13... with the maximal number of non-zero charge densities Qj j 1 ...j p obeying (4.15). 

5. Cosmological solutions to the field equations for D=5, 9, 13, ... 

Here we give explicit examples of cosmological type solutions for the dimensions from (4.22) 
when the non-distinguished coordinates are all space-like, i.e. 

£l = ... = e n = l. (5.1) 
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First we show that all scale factors are the same up to constants: 



ft{u) = ftu) + c\ (5.2) 

From the definition of running constants (4.2) and the (anti-) self-duality of the charge density 
form (4.15), it follows that 

Y<P = Y& + const, (5.3) 
iei jei 

where I is an arbitrary subset of Iq = {1, . . . , n} of length n/2 = 2m and 

I = I \I, (5.4) 

is "dual" set. For D = 5 case eqs. (5.3) read: (j) 1 + 4> 2 = 4> 3 + (j) A + const and two other 
relations obtained by permutations. 

From relations (5.3) one can see that all ft should coincide up to constants (for i £ I 
and j € / it is sufficient to consider another equation with I\ = (I \ {i}) U {j}) instead of I 
in (5.3) and find from both equations that ft and ft coincide up to a constant). 

Thus, we are led to (5.2). In what follows we put c l = which may always be done via 
a proper rescaling of y-coordinates. This also implies that non-running charge density form 
Qi ...i p is self-dual or anti-self-dual in a flat Euclidean space R n , i.e. 

Qi ...i p = i ^ _|_ ^ , £ i . . .i p jo ■ ■ -j p Q J ° ' ' ' Jp = ±(*Q)i ...ip- (5-5) 

The Lagrangian and total energy constraint are given by 

Lq = \G t] ftft -V Q + ^\ (5.6) 
E Q = l -G l3 ftft + V Q + hp 2 = 0, (5.7) 



with the potential being 



V Q = \Y1 exp ( 2 Y, ~ 2X A > ( 5 - 8 ) 

/ V kei J 



see (3.13)-(3.15). 

The field equations for <f> and ip from Lq are 



Y ^ft + Y Q 2 ^ ex p 2 E ^ - 2 v = °> ( 5 - 9 ) 

j=i I V fee/ / 

(^ + ^Q 2 (/)(-A)exp ( 2Y^ k -2Ayj J =0. (5.10) 
/ V fee/ / 
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Since all <p l satisfy (j? = <\> one finds YTj=\ = X)?=i(^»j ~~ = 0- ~ n ) ( t ) - Next, denning 



and noting that 



Q 2 ^Q 2 (/)/0 



for any i = 1, ...,n, one finds that the field equations (5.9) and (5.10) become 

Q 2 exp(n0 — 2Xip), 



r 2(n-l) 
i/3 = AQ 2 exp(n0 — 2A99) 



(5.11) 
(5.12) 

(5.13) 
(5.14) 



Finally, since the exponents in (5.13) and (5.14) are the same these two equations can be 
combined into one as 



with the definitions 



n 



f = M 

A = —K, K = \ 2 

2 4(n - 1) 



n 



assumed. 

The first integral of (5.15) is 



where C is an integration constant. 
Let K / 0, or 



A 2 # 



n 



4(n - 1) 



- \ 2 



Equation (5.15) has several solutions: 

/ = -ln \z\2A\ 1 ! 2 



with 



z = —= sinh 

Vc 
1 

: sin 



v-c 

U - Uq, 



(u - u )VC 
(u - Uq)^C 



A < 0, C > 0; 



A < 0, C < 0; 
A < 0, C = 0; 
= -j= cosh (it - i^v 7 ^] , A > 0, C > 0. 



(5.15) 

(5.16) 
(5.17) 



(5.18) 

(5.19) 
(5.20) 

(5.21) 

(5.22) 
(5.23) 
(5.24) 
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One can relate the solutions given in (5.20) and (5.21) - (5.24) to (f) and (p by using (5.13) 
to construct the following relationship 



<p = 2(n - l)\<f>, (5.25) 

which has the solution 

<p = 2(n - l)\<f> + C 2 u + d, (5.26) 
where C 2 ,Ci are integration constants. Combining (5.26) with (5.16) gives 

1 [A(C 2 « + Ci) + /(«)], <p= - n — (C 2M + Ci)-^^. (5.27) 



r 2(l-n)K L v r 4(1 -n)^T v ±y if 

Applying this to the zero energy constraint 



and using (5.18) we get 



or equivalently, 



E Q = ±n(l - n)tf + \p 2 + IqV/M = (5.28) 



C n(C 2 ) 2 



K AK{n - 1) 



n 



This tells us that only the three cases (5.21), (5.23) and (5.24) occur when real solutions 
are considered. The solution (5.22) will be considered in a possible future work with a pure 
imaginary scalar field and A (this is equivalent to a "phantom" field that may support the 
so-called bouncing solution). 

Collecting these results together the solutions for the metric, scalar field and (p + 2)-form 

are 

n 

ds 2 = we 2n ^du 2 + e 2 ^ J^(*/) 2 (5.31) 

i=l 

F = e 2 ^duAQ, Q = _L^Q i0 . ^dy^ A . . . A (5.33) 

with <j){u) given by (5.27) and the function f(u) given by (5.20) and (5.21), (5.23), (5.24). 
Here the charge density form Q of rank n/2 = 2m is self-dual or anti-self-dual in a flat 
Euclidean space M n : Q = ± * Q, the parameters C 2 , C obey (5.30) and the dilatonic coupling 
constant A is non-exceptional, see (5.19). 
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5.1 Special attractor solution for C = 

Here we examine some properties of the simplest cosmological type solution given in (5.23). 
For this solution A < and hence 

A 2 < Aq. (5.34) 

Since C = for (5.23) one has C2 = from the condition (5.30). Finally, without loss of 
generality we put uq = 0. 

To get a physical understanding of this solution one should change the "time" coordinate 
u to the proper time coordinate r. In order to get the correct sign for the proper time we fix 
the sign in the relationship between u and the proper time as 

d T = -e n(f,(u) du, (5.35) 

where 

<P = 2 (n-l)K Hul2All/2) > A = A ^~ 2XCl - (5-36) 
Integrating (5.35) and taking a suitable choice of reference point we get 

\a\\2A\ l ' 2 T = {u\2A\ 1 / 2 ) a , (5.37) 

where u > and 



a = < °- ( 5 - 38 ) 



A 2 + A| 

Since a < 0, r = t(u) is monotonically decreasing from infinity when u = to zero when 
u = 00. 

The metric (5.31) now reads 

n 

ds 1 = wdr 2 + Bt 2u ^(d^) 2 , (5.39) 
i=i 

where r > and 

■>= n + J(n-iy B =<HM 1/2 )" <"0) 

By putting A = and w = — 1 in the above solution we get a cosmological powerdaw 
expansion with a power v = 2/n that is the same as in the case of D = 1 + n dust matter with 
a zero pressure (see for example [18, 19, 16]). This is not surprising since it can be argued 
that the collection of branes with charge densities obeying (anti)-self-duality condition (5.5) 
behaves as a dust matter. Indeed, we know that in the absence of a scalar field the solution 
with a single brane described by a set I C {1, n} is equivalent to an anisotropic fluid with 
equations of state pi = — p for i G / and Pi = p otherwise, i = 1, . . . ,n [5]. Here pi is a 
pressure in the i-th direction and p is the energy density. Thus, our collection of branes is 
equivalent to a multicomponent fluid [20] with coinciding densities, since all Q 2 (I) are equal 
and <p % = (p. For any i the collection of branes can be split into pairs with sets /, I such 
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that i G I. The first brane gives pressure pi = —p and the second one gives p\ = p (energy 
densities are the same ). Hence the total pressure is zero for the pair and, thus, for the whole 
collection of branes. 

Finally, we note that the solution (5.39) is attractor solution in the limit r — ► +00, 
or u — > +0, for the solutions with A < given by (5.21). This follows just from relation 
sinhn ~ u for small u. 



5.2 Kasner-like behavior for r — ► +0 

We now consider u — ► +00 asymptotical behavior of solutions with i) sinh- and ii) cosh- 
functions corresponding to (5.21) and (5.24), respectively. In both cases we find that the 
asymptotic behavior is Kasner like i.e. the metric and scalar field take the form 

n 

ds 2 as = wdT 2 + Y,T 2al A i {dy i ) 2 , ip as = a ip \n T + ip , (5.41) 
i=i 

where Ai > 0, ipo are constants. The Kasner parameters obey 

n n 

= J>*) 2 + K) 2 = 1. (5.42) 
i=i i=i 

In the first case i) with sinh-dependence and A < the proper time r is decreasing when 
u — > +00. From eqs. (5.20), (5.21), (5.24), (5.27) and (5.35) we find the following asymptotic 
behavior as u — ► +00 

<f> ~ -bu + const, 6 = — — > 0, (5.43) 

2(n — 1)K 

ip ~ (-A^C 2 + A^C)^^ + const, (5.44) 
t ~ const exp( — nbu). (5.45) 

Using these asymptotic relations and writing everything in terms of proper time one find that 
the metric and scalar field take the following asymptotic forms 

n 

ds 2 as = wdT 2 + r^ n A Y J (.dy 1 ) 2 , tpas = a^lnr + <po, (5.46) 



aj = l--, (5.47) 



as r — > +0. Here > 0, 990 are constants. The relationship for q 2 comes from (5.44), 
(5.45) and agrees with (5.42) and a 1 = l/n. Using (5.43)-(5.45) one can obtain the following 
relationship: signfa^] = — signfCy. 

In the second case ii) with cosh-dependence and A > the proper time r decreases as 
u — > +00 for AC2 > and increases for AC2 < 0. In this case we also get an asymptotical 
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Kasner type relations (5.46)-(5.47) in the limit r — > +0 with sign[a^] = — sign [A]. In both 
sinh— and cosh — cases the Kasner sets a = (a 1 = 1/n, a v ) obey the inequalities 



for all brane sets I. This is in agreement with a general prescription of the billiard represen- 
tation from [13]. In the case ii) we obtain the Kasner type behavior (5.46)-(5.47) in the limit 
r — ► +00 with signfa^] = sign [A]. In this case 



for all /. In case ii) the scalar field has a bouncing behavior in the interval r £ (0, +00). 
6. Conclusions 

In this article we have examined a system where D = n + 1 dimensional gravity was coupled to 
a scalar field plus ap + 2 rank form field. The ansatz employed here was to take the metric as 
diagonal, and the rank p + 2 form field to have a composite electric Sp-brane form. All ansatz 
functions depended only on the one distinguished coordinate, u. Under these conditions the 
initial model could be reduced to an effective 1-dimensional cr-model which greatly simplified 
the study of the system. The diagonal character of our metric ansatz resulted in there being 
constraint equations among the charge densities of branes associated with the p + 2 form 
field. By examining these constraint equations we showed that for certain odd values of the 
spacetime dimension, given by D = Am + 1 = 5,9,13..., the system allowed the maximal 
number of the charged branes. We also proved the absence of such maximal configurations 
for p = 1 and even D. 

For these special odd dimensions and non-exceptional dilatonic coupling (A 2 7^ 477^1) ) 
we wrote down exact solutions given in equations (5.31)-(5.33). We examined the simplest of 
these solutions given by (5.23). On converting the distinguished coordinate u to the proper 
time r, we found this solution corresponded to a cosmological power-law expansion. By 
taking the limit of vanishing dilatonic coupling, A = 0, we showed that this solution reduced 
to a cosmological model with dust matter. On physical grounds this is exactly what one 
would expect in this limit. We also investigated an asymptotical Kasner type behavior of the 
solutions for small (r — > +0) and large (r — > +00) values of proper time. 

Directions for possible future work include: (i) the investigation of the properties of the 
solution in (5.22) using pure imaginary scalar field and A (with a hope to investigate the 
possibility of bouncing behavior); (ii) the investigation of static, non-cosmological solutions, 
when w = +1 and the distinguished coordinate u is spacelike; (iii) examining the exceptional 
case, when K = (i.e. A 2 = 4 ( n n _ 1 ) ) in (5.17). 

Using the results obtained in this paper and in [13] we will present in future work [29] 
some examples of breaking of the never-ending oscillating behavior near the singularity by 




(5.48) 



f/ 7 (a) < 



(5.49) 
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constraints. Physically this may indicate that the form of the ansatz considered here (diagonal 
metric, composite electric ansatz for the form field, etc.) may prevent the occurrence of the 
"chaotic", never-ending oscillating behavior (see [13, 14, 15, 21, 22] and references, therein). 

It should be noted that recently N.S. Deger has obtained in [23] some 5-brane solu- 
tions in D = 11 supergravity with non-standard intersection rules. These solutions contain 
"dangerous" intersections of branes and give certain examples of solutions with charge den- 
sities obeying constraint equations. It seems that the "non-standard" solutions from [23] 
with vanishing Chern-Simons term may be obtained as a special case of the so-called "block- 
orthogonal" solutions [24, 25, 26]. In these block-orthogonal solutions the branes split into 
several blocks. The branes belonging to one block behave as one "effective" brane. The 
branes belonging to different blocks obey the "standard" intersection rules [12, 27, 28]. In 
future work [29] we will verify whether the solution from our paper may be obtained using 
the modified one-block solution from [25, 26]. 
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